Gamma-rays Prom Warm WIMP Dark Matter Annihilation 
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The weakly interacting massive particle (WIMP) often serves as a candidate for the cold dark 
matter, however when produced non-thermally it could behave like warm dark matter. In this paper 
we study the properties of the 7-ray emission from annihilation of WIMP dark matter in the halo 
of our own Milky- Way Galaxy with high resolution A''-body simulations of a Milky- Way like dark 
matter halo, assuming different nature of WIMPs. Due to the large free-streaming length in the 
scenario of warm WIMPs, the substructure contend of the dark matter halo is significantly different 
from that of the cold WIMP counterpart, resulting in distinct predictions of the 7-ray signals from 
the dark matter annihilation. We illustrate these by comparing the predicted 7-ray signals from the 
warm WIMP annihilation to that of cold WIMPs. Pronounced differences from the subhalo skymap 
and statistical properties between two WIMP models are demonstrated. Due to the potentially 
enhanced cross section of the non-thermal production mechanism in warm WIMP scenario, the 
Galactic center might be prior for the indirect detection of warm WIMPs to dwarf galaxies, which 
might be different from the cold dark matter scenario. As a specific example we consider the 
non-thermally produced neutralino of supersymmetric model and discuss the detectability of warm 
WIMPs with Fermi 7-ray telescope. 

PACS numbers: 95.35.-|-d,95.85.Pw 



I. INTRODUCTION 

The so-called dark matter (DM) , discovered ~ 80 years 
ago in the astronomical observations, is still one of the 
biggest mysteries in the fields of physics, astronomy and 
cosmology. To understand the nature of DM particles 
is a big challenge of the community. There are several 
ways being proposed to detect the WIMP DM particles 
(see e.g., [l|), among which the indirect search through 
the cosmic ray (CR) particles is the most active one in 
recent years due to the operation of several new gener- 
ation satellites, such as PAMELA, Fermi and AMS02. 
In many kinds of CR particles, the anti-particles, 7-rays 
and neutrinos are good probes to search for DM signals. 
Especially, 7-rays are widely discussed, due to the simple 
propagation and the high sensitivity detections from both 
spatial and ground-based telescopes. The constraints on 
the DM parameters become stronger and stronger in re- 
cent years thanks to the Fermi 7-ray observations [2t6| . 

One of the key problems in the study of the 7-ray emis- 
sion from the WIMP DM annihilation is the density dis- 
tribution of DM. It is observationally very difficult to 
determine the density distribution of DM, especailly at 
small scales. Currently the postulated best knowledge 
about the DM density distribution comes from the nu- 
merical N-body simulations (e.g., [3-Q)- 

The initial matter power spectrum which describes cos- 
mic density perturbation depends on the particle nature 
of DM. For the cold DM (CDM), the particle velocity 
when decoupling is negligible and the corresponding free- 
streaming length is very short. The small free-streaming 



length enables structures down to very small scales to 
form. 

The CDM scenario has been shown to be in good agree- 
ment with the observations of the cosmological large scale 
structures. However, it has been a long time problem of 
the CDM scenario that the expected structures are in- 
consistent with observations at sub-galactic scale (e.g., 
[I(]|4l4l | ) . One possible solution of this problem is the 
warm DM (WDM) scenario ( |15l - [l7| . or a recent review 
[18|). In general, with a thermal distribution, the parti- 
cle mass of the WDM should be as light as '-^keV. After 
decoupling the velocity of WDM can be fast enough to 
introduce a large free-streaming scale below which the 
structures are smoothed out. Thus the formation of small 
scale structures in the WDM scenario can be suppressed. 

If the DM is finally proven to be warm, the impact on 
the detection of DM particles is fatal, because most of 
these experiments aim to search for the weakly interact- 
ing massive particles (WIMPs) which are traditionally 
cold. For the canonical WIMPs, when produced ther- 
mally in the early Universe, the velocity is non-relativistic 
after decoupling and they behave like CDM. Alterna- 
tively the WIMPs if produced non-thermally, it can be 
warm [19| - |22| |. In Ref.[20] the authors showed explicitly 
that the power spectrum of this non-thermally produced 
WIMPs has a clear suppression at small scales. The non- 
thermally produced WIMP scenario will have some inter- 
esting properties for the indirect .search of DM, because 
1) compared with the light WDM the mass of the non- 
thermal WIMPs lies within the range of most high en- 
ergy CR detectors, and 2) in contrast with the thermally 



produced WIMPs, the annihilation cross section of non- 
thermal WIMPs can be larger due to the lack of direct 
constraints from the relic density. We will discuss the 
possible 7-ray signatures from such non-thermal WIMP 
DM annihilation in this paper. 

In this paper we focus on predicted DM annihilation 
signals from the Milky Way halo and its substructures 
based upon high resolution simulations of WDM in [2j] . 

This paper is organized as follows. In Sec. II we 
briefly introduce the picture of the non-thermally pro- 
duced warm WIMPs. In Sec. Ill we describe the numer- 
ical simulations used in this work and the DM density 
distributions for the smooth halo and subhalos according 
to the simulations. The signatures of 7-ray signals and 
detectability analysis are discussed in Sec. IV. Finally 
Sec. V is the conclusion. 



II. NON-THERMALLY PRODUCED WARM 
WIMPS 

The DM particles can be non-thermally produced by 
the decm^ of topological defects such as cosmic string 
p^A I20I |24| . |25| . For example, we consider a model with 
an extra U(l) gauge symmetry which is broken by the 
vacuum expectation value 77 of a scalar field S [l9| . Cos- 
mic strings will be formed during the symmetry break- 
ing phase transition taking place at the temperature of 
Tc ^ rj. After the transition, the infinite long string net- 
work coarsens, and more closed string loops form from 
the reconnection of the long strings. The tension of the 
cosmic string is determined hy ^ ^ rf. Cosmic string 
loops lose their energy dominantly through gravitational 
radiation. When the radius of a loop becomes to the or- 
der of the string width, the loop will self-annihilate into 
its constituent field, such as scalar boson S. The DM 
particle x can be produced by the decay of these heavy 
particles. 

When the temperature of DM is higher than the freeze- 
out temperature T^{^ O(GeV)), DM particles produced 
by cosmic string loops still keep in chemical equilibrium 
with standard model particles. Only the DM particles 
produced below T^ will contribute to the non-thermal 
DM relic density JInt- It is found that the non-thermal 
DM is mostly contributed by the loops decaying at T^ 
(Eq. (A2)). Therefore the DM production process does 
not affect the big bang nucleosynthesis (BBN) results. 
Through adjusting the model parameters, the relic den- 
sity of DM can also be naturally explained (for more 
details, see Appendix A). 

In such scenario, DM particle x may carry large initial 
momentum pc due to the decay of the heavy particle. 
Pc can be written as pc — aTc where a is a numerical 
factor determined by detailed model. Here we define a 
typical variable re = a(t)p(t) /m-^ which is a constant 
during the cosmic evolution [20[ . If we choose the cosmic 
scale factor at the present time a{to) = 1, Tc can be 
understood as today's velocity of the DM particles if no 



structure formation. The comoving free-streaming scale 
Rf is given by ^ 






2r-ctEQ(l + 2;EQ)^ln 
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(1 + ^eq)2 
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^c(1 + ^eq). 
where "EQ" denotes the radiation-matter equality. 
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FIG. 1: Linear matter power spectra of CDM (red, short- 
dashed), canonical light WDM (blue, long-dashed) and non- 
thermal warm WIMP (black, solid). 

The free-streaming of DM particles will imprint on the 
late time structure formation. This effect can be sim- 
ply seen by the matter power agectrum of DM. We use 
a modified version of CAMB^ [2a| to calculate the mat- 
ter power spectrum of the non-thermally produced DM 
scenario, shown in Fig. [T] Here we adopt re = 10^^. 
Note the mass of the non-thermal DM does not explic- 
itly affect the calculation of the power spectrum because 
its effect can be cancelled by the initial momentum (see 
the definition of re). For comparison the power spectra 
for CDM and the canonical light WDM are also shown. 
The power spectrum of the canonical WDM corresponds 
to a sterile neutrino with mass ~ 2 keV, which is also the 
input power spectrum of the N-body simulation (see be- 
low Sec. HI). We can see that a clear suppression of the 
power at small scales appears both for the light WDM 
and the warm WIMP scenarios. The free-streaming prop- 
erty makes the non-thermal DM behave similarly with 
WDM. Due to the similarity of the input power spectra 
of the light WDM and non-thermal warm WIMPs, we use 
the simulation results for the light WDM in the following 
discussion of the indirect detection of warm WIMPs. 



^ http://cairLb.info 



III. NUMERICAL SIMULATION RESULTS 

In this section, we describe briefly numerical simula- 
tions used in this work and present the properties of 
the DM distribution based on the numerical simulations. 
The simulations used in this study are two matched ultra- 
high resolution simulations of a Milky sized DM halos run 
with different nature of DM models but with otherwise 
same numerical setup as well as cosmological parame- 
ters. For the CDM simulation, we use "Aq-A-2", from 
the Aquarius Project [7|. In order to facilitate compar- 
ison of DM annihilation emission from cold and warm 
DM models, for the same halo, we further performed a 
high resolution simulation assuming WDM model by us- 
ing the same phase in the initial density field as that of 
the "Aq-A-2" simulation but a different matter spectrum 
matching a particular WDM model. In the numerical 
calculation of this paper, we adopt a 2 kev sterile neu- 
trino [27| as our WDM model which lies within bound of 
Lya constraint [27j . The chosen WDM introduce a cutoff 
emerging at a wavenumber k ~ lOh Mpc~^ in the initial 
matter power spectrum, below which the power spectrum 
is well consistent with that of CDM [2^ . In the scenario 
of non-thermal WIMPs, such as a heavy particle S de- 
caying into two WIMPs Xj for x around 100 GeV, it 
requires the mass of particle S around 10* GeV [20|. In 
our simulation, the mass of the "particle" is 1.37 x lO** 
M0, and the number of "particles" is larger than 100 mil- 
lion within r2oo, the radius inside which the mean DM 
density is 200 times of the critical density. Therefore the 
lowest mass subhalos resolved in our simulation is 3 x 10^ 
M0 if requiring more than 20 "particles" for a subhalo. 
The total mass within r2onof the halo is about 1.8 x 10^^ 
M0. See Table 1 of Ref. [2j] for the basic information of 
the simulations. For a more detailed description of our 
simulation, please refer to Ref. [23] . 



A. Smooth halo 

The density profile of the smooth component of the 
simulated halo of the CDM simulation was analyzed in 
Ref. 125. It was shown that the smooth halo density 
profile can be well fitted with an Einasto profile [22 



p{r) = p_2 exp 



r-2 



-1 



(2) 



where p_2 « 0.14 GeV cm ^, 7'_2 ~ 15.7 kpc and 
a « 0.17 [Hi. The local density of DM is then given 
as pq « 0.44 GeV cm~'^ at Rq = 8.5 kpc. A higher local 
density compared with the canonical 0.3 GeV cm~'^ was 
also found in recent studies |30l - l32 |. 

For the WDM halo the density profile of the smooth 
halo is essentially the same as that of CDM down to the 
numerical resolution limit of our simulation [33|. The 
expectation that a core may appear in the center of the 
halo for WDM due to phase space density constraint [3J- 
|36| is not clearly seen in the simulation, this is because 



the core size of the Milky Way sized halo is predicted 
to be smaller than resolution limit of our simulation and 
thus is not resolved. It was shown recently that a density 
core was indeed observed in WDM simulations, at a scale 
smaller than 100 pc for 1 — 2 keV WDM and halo mass 
10* - IQio Mq [33. For the Milky Way like halo the 
expected core will be even smaller, and the halo density 
profile will be indistinguishable from that of CDM halo, 
within the precision of sub-degree of the present 7-ray 
detectors. In this work we adopt the same equation ^ to 
describe the density profile of the smooth halo for WDM. 



B. Subhalos 

Based on the simulation results, we find 20529 gravi- 
tational bounded subhalos for CDM simulation and 219 
subhalos for WDM^ simulation within the virial radius 
of the main halo. The minimum mass of the resolved 
subhalo is found to be ~ 3 x 10^ Mq, and the maximum 
mass is about 10^° M©. 

We define the annihilation luminosity of a subhalo as 
Li = J pfdVi. In the work we adopt the Navarro- Frenk- 
White (NFW, [H]) profile for the subhalos. The determi- 
nation of the parameters of the NFW density profile from 
the simulated circular velocity profile can be found in the 
Appendix of Ref. [33| • For WDM subhalos we employ a 
constant density core with size Vc ~ 0.03 x ( j^^'^-i ; 
kpc, where a is the velocity dispersion of the subhalo [4C 



-0.5 



Beyond re the density distribution is identical with NFW 
profile. The 7-ray flux from DM annihilation of this sub- 
halo is then proportional to Li/df, with di the distance 
of the subhalo from Earth. To calculate di of each sub- 
halo, a random location of the solar system which is 8.5 
kpc away from the halo center is chosen. 

The mass-luminosity and mass-fiux scattering plots of 
the subhalos are shown in Fig. [2] From the mass-flux 
relation we see that in general, subhalos in CDM case are 
more brighter than that of WDM because of relatively 
lower concentration of subhalo in WDM comparing to 
CDM Hi]. There are also fewer subhalos of WDM which 
can have comparable fluxes as that of CDM. Especially 
we find the most massive subhalos are usually not the 
brightest objects. The subhalos with masses 10^ — 10^ 
Mq have larger probability to give high fluxes |41| . 

For the CDM case it is expected that there should be a 
large number of unresolved substructures below the res- 
olution limit of the simulation, which can extend to a 
mass comparable or even lower than that of the Earth, 
lO"*" M© m, HI]. To include the contribution of unre- 
solved subhalos we have to extrapolate the subhalos to 
lower mass, according to the statistical properties of the 
resolved subhalo distributions [7|, |4J] . We present the ba- 



^ Note for WDM case, the number of subhalos might be over- 
estimated due to the numerical fragmentation of filaments 13^ . 
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FIG. 2: Annihilation luminosity (L, left) and relative flux (F, right) versus mass of subhalos for CDM and WDM simulations. 



sic statistical results of the subhalos of CDM and WDM 
based on the simulations in the Appendix. For the WDM 
case, because free streaming length of the chosen WDM 
particle is as large as 200 kpc, the smallest dark mat- 
ter halo expected to form in the model is therefore about 
2.5 X 10^ M© [13], corresponding to ~ 10^ particles in our 
simulation and hence are well resolved in our simulation. 
Thus we believe that our WDM simulation has resolved 
all subhalos and theorefore has no unresolved subhalo 
compounent. There are some spurious subhalo formed 
in our simulation via artificial fragmentation of filaments 
as noted by Ref. [33|, however we expect that contribu- 
tion to the annihilation luminosity due to these spurious 
subhalos is small because of their low abundance. We do 
not consider them in the following analysis. 



C. J-factors 

The 7-ray fiux observed at the Earth from DM anni- 
hilation can be written as 



(j)^{E^,il}) 



pIRq (av) dN 



47r 



2m2 dE-, 



X Ji^P), 



(3) 



where m^^ is the mass of the DM particle, (av) is the 
annihilation cross section weighted with the velocity of 
DM particle, 4^ is the 7-ray yield spectrum per annihi- 
lation. The dimensionlcss astrophysical J-factor, related 
to the DM density profile, is defined as 
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where ip is defined as the angle between the observational 
direction and the Galactic center direction for observer 
at the Earth. The integral is done along the line of sight 
(LOS). Taking the detectors angular resolution into ac- 
count the J factor for a resolved subhalo is defined as 



JLhW 
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where Li , di and ipi are the luminosity, distance and cen- 
tral direction of the ith halo. The exponential term on 
the right hand side corresponds to a Gaussian smooth 
with width a. 

Based on the numerical simulation of WDM wc cal- 
culate the J factor of the smooth halo and the subhalos. 
The skymaps of the J factors of the smooth halo, resolved 
subhalos and the total result for WDM are shown in Fig. 
|3l The colorbar shows the value of log( J). For resolved 
subhalos we employ Gaussian smoothing with a = 0.5°. 

The skymaps of the CDM subhalos based on the sim- 
ulation Aquarius has been given in Ref. [8|. To com- 
pare with the skymap of WDM subhalos given in this 
work we have also shown the skymaps of CDM in Ap- 
pendix (see Fig. [TU)) . From those two figures we can 
clearly see the differences between the CDM and WDM 
annihilation signals from subhalos. For CDM, there is 
a non-negligible diffuse component from the unresolved 
subhalos, especially at the directions far away from the 
Galactic center. The number of the potentially visible 
subhalos above the diffuse component is much higher for 
CDM than WDM. There is also difference in the Galactic 
center due to expected presence of a core in the WDM 
scenario. However, we may over-estimate the size of the 
core in this work compared with that found in the simu- 
lations [33] ■ The actual difference may be smaller. 

The accumulative subhalo number which represents 
the subhalos with J factor greater than some value Jsub 
versus Jgub is shown in Fig. Ul The different lines in 
each group represent the random choice of the location 
of the solar system in the halo, with distance fixed to be 
8.5 kpc from the center. We can see that the number 
distribution of WDM is flatter than that of CDM. This 
is because in the CDM case the relative weight of smaller 
subhalos compared with larger ones are higher than that 
in the WDM case. The property presented in Fig. SI 
if detectable, is useful to probe the nature of the DM 
particles. 




FIG. 3: Skymaps of the J-factors of the main halo (top-left), resolved subhalos (top-right) and total contribution (bottom) for 
WDM. 




FIG. 4: Accumulative number versus J of subhalos. 



IV. GAMMA-RAY SIGNALS 

In this section wc study the 7-ray signals from the 
warm WIMP annihilation. We will present the astro- 
physical 7-ray background and the detectability of the 
7-rays from warm WIMP annihilation by Fermi. 



A. Benchmark models of supersymmetric DM 



ing the constraints from e.g., 7-rays and antiprotons, the 
cross section can not be arbitrarily large. The constraint 
from PAMELA antiproton data showed that the allowed 
boost factor'^ of neutralino-like DM should be less than 
10 for O(IOO) GeV DM ^. The new constraints from 
Fermi observations of dwarf galaxies also gave an allowed 
boost factor of several for O(IOO) GeV DM H]. Taking 
the above contraints on the WIMP annihilation cross sec- 
tion into account, we give two explicit benchmark models 
to realize the cold and warm WIMP scenarios in super- 
symmetric DM models. 

In the supersymmetric (SUSY) theory with R-parity 
conservation, the lightest neutralino, which is the combi- 
nation of gaugino and higgsino, is a well-motivated can- 
didate of DM |4g. In general, there are four param- 
eter regions to obtain the correct thermal relic density 
of neutralino: (1) all the sfermions are light, neutralinos 
annihilate via t-channel sfermions exchange; (2) neutrali- 
nos scatter with sfermions with nearly mass degeneracy 
which is so-called "co-annihilation"; (3) Xi has signifi- 
cant component of Higgsino or wino, with the main an- 
nihilation channel to heavy gauge boson or Higgs; (4) 
neutralinos annihilate via s-channel Higgs resonance with 
2to^o ~ 771^0 , or nifio , mjjo . In the first region, the light 
sfermions are stringently constrained by recent LHC re- 
sults [43, 13 • I'^ tti^ "co-annihilation" region, the neu- 



For the warm WIMP, the annihilation cross section 
may be larger than that of cold WIMP which are con- 
strained by the relic density of DM. However, consider- 



3 Defined as ((tii)/3 x IQ-^e 



TABLE I: Relevant parameters for the two benchmark models. The unit of mo, m„ , m„ , mi/2, Aq, m^o is GeV, and of 
(av) is cni"'s~"'". 





mo 


"^H^ 


m„ 


»Tll/2 


Ao tan/3 sign(/^) 


m-a 


{av) 


Warm WIMP 


1200 


1300 


788 


500 


-1000 40 + 


211 


2.70 X 10^^^^ 


Cold WIMP 


1200 


1300 


824 


500 


-1000 40 + 


211 


1.38 X lO"^'' 



tralino annihilation cross section is often much smaUer 
than the "natural value" 3 x 10^^^ cm^s^^ Thus it is 
difficult to observe the products of DM annihilation in 
indirect detections. In the third region, neutralino an- 
nihilation could produce large flux of 7-rays due to cas- 
cade decay of gauge boson or Higgs. However, a signif- 
icant component of Higgsino or wino in the neutralino 
might induce large interaction between DM and nucleon, 
which is stringently constrained by recent direct detec- 
tions, such as XENONIOO [H. 

Here we consider two benchmark models in the "Higgs 
funnel" region as the cold and warm WIMP candidates. 
The DM annihilation is enhanced by s-channel pseu- 
doscalar Higgs exchange with resonance effect tti^o ^ 
2m^a , and the main final states of DM annihilations are 
66". The ATLAS and CMS collaborations have discov- 
ered a 125 GeV Higgs-like boson [50, [5l|- Because the 
Higgs in MSSM is lighter than Z boson at the tree level, it 
requires large stop mass parameter or large mixing term 
to acquire corrected Higgs mass. It can be interpreted 
by some particular parameter configurations. Since here 
we employ the benchmark models as illustration and em- 
phasize the difference between warm and cold WIMPs, 
we do not consider this issue of Higgs mass in this work. 

To acquire a moderate A^ mass easily, we consider the 
"NUHM" scenario [52, llj] , in which the Higgs mass pa- 
rameters niH^ j^^ffd £^t GUT scale are different from other 
scalar masses toq- The particle spectrum, DM thermal 
relic density and annihilation cross section for the bench- 
marks models are calculated by SuSpect fSJ] and mi- 
crOMEGAs [5^, l5y| , and are summarized in Table HI For 
the warm WIMP model adopted here, the thermal relic 
density of DM is f^/i^j^ ~ 4.33 x lO^^^, much smaller than 
the observational value ^h? ^ 0.11. Therefore it must 
be produced via non-thermal mechanism (sec Appendix 
A and [241 )■ Given the particle models of DM, the 7-ray 
spectrum from the decay and fragmentation of the final 
state particles is calculated using the PYTHIA simula- 
tion tool IstI. 



In the benchmark 
i^ « 211 GeV, {av) 



models the DM parameters are 
w 1.38 x 10-26 cm^ s-i for cold 



* The potential phenomenology problem of this region may bo the 
large contribution to rare meson decay, such as B^ — > Xs'^ and 
Bs — )■ fJ-'^ tJ.~ , due to light pseudoscalar Higgs and large fan/3. To 
avoid violating meson decay observations, some special parame- 
ters are needed to suppress total SUSY contributions from Higss 
sector and chargino-squark sector. 



WIMP and 2.70 x 10-^5 cm^s"! for warm WIMP, and 
the annihilation final state is about 86% 66 -f IA%t^t~. 
The cross section for warm WIMP corresponds to a boost 
factor of 9, which is roughly compatible with the present 
constraints from indirect detections. Note these con- 
straints are applicable for neutralino DM. For other an- 
nihilation final states such as the leptons the constraints 
might be different, and even larger boost factor could be 
possible. 



B. Astrophysical background 

To discuss the detectability of DM we have to take the 
astrophysical background of diffuse 7-rays into account. 
We use the GALPROP^ Q code to calculate the Galac- 
tic diffuse 7-ray background. The propagation parame- 
ters adopted are: Dq = 6.59 x 10^* cm2s-\ S = 0.30, 
VA = 39.2 km s~^. Zh = 3.9 kpc, according to the fit 
to the B/C data [53 ■ The injection spectra of nuclei are 
adopted as 7" = I.9I/72 = 2.40 for rigidity below/above 
10 GV, which can basically reproduce the recent mea- 
surements of proton and Helium spectra by PAMELA 
[601, as shown in [61|. Note, however, this simple injec- 
tion model may not well describe the detailed harden- 
ing structures of the CR spectra above several hundred 
GV, or the difference between proton and Helium spec- 
tra ( [60, [63, ll^ ) ■ For CR electrons, the injection spec- 
tra are ^f = I.5O/7I = 2.56 for rigidity below/above 
3.55 GV as derived according to the pure background fit 
to the newest e^e^ data [6l|. Such a pure background 
component cannot explain the e~^e~ excesses revealed by 
several experiments p34| - [68| . As illustrated in [601 ^^^ 
contribution to the total diffuse 7-rays from the extra 
astrophysical sources of e+e", e.g., pulsars, is always 
negligible in all regions of the sky. For the purpose of 
the current study, we think it is enough to employ such 
a rough background model. Finally the extra-galactic 
7-ray background is adopted to be the Fermi measured 
results, ^G ~ 5.89 x 10-^(£;/GeV)-2''4 cm "^ s'^ sr"! 
GeV-i jifl. 

We calculate the total 7-ray skymaps for the cold and 
warm WIMP scenarios based on the two SUSY bench- 
mark models given the previous subsection. The total 
7-ray skymaps above 10 GeV of both the astrophysical 



http://galprop.stanford.edu/ 




FIG. 5: Skymaps of the total 7-ray emission with background predicted by GALPROP included, for energies E > 10 GeV. 
The left panel is for the cold WIMP case, and the right panel is for the warm WIMP case. 



background and the DM contribution arc shown in Fig. 
[5l The left panel is for cold WIMP and the right panel 
is for warm WIMP respectively. The detectability of the 
DM signal in presence of the astrophysical background 
will be discussed in the followings two subsections. 



Gamma-rays from warm WIMP annihilation: 
Galactic center 
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FIG. 6: Gamma-ray spectra in 20° x 20° region around the 
Galactic center for cold and warm WIMP scenarios. Shaded 
region represents the expected background of the GALPROP 
model (see the text). 

Fig. [6] shows the expected 7-ray spectra of the WIMP 
annihilation in 20° x 20° region around the Galactic cen- 
ter. For comparison we also plot the calculated diffuse 
background described in Sec. IV. B. There are all-sky 
survey data of diffuse 7-rays from Fermi, available from 
the Fermi Science Support Center^. It was shown that 
the GALPROP model could reproduce the observational 



data within a precision of factor 2 [7l|. Therefore here 
we simply employ the model results for comparison. An 
uncertainty of 2 times of the GALPROP background is 
represented by the shaded region. 

It is shown that for warm WIMP scenario we may ex- 
pect larger flux of 7-rays in the Galactic center, sim- 
ply due to a larger annihilation cross section of warm 
WIMPs. Up to now there is no clear evidence to show the 
existence of signals from DM in the Fermi data^. How- 
ever, we may expect that the warm WIMP scenario could 
have better detection perspective than the cold WIMP 
scenario. 



D. Gamma-rays from warm WIMP annihilation: 
subhalos 



Finally we discuss the detectability of DM subhalos. 
Fig. [7] shows the integral fluxes above 100 MeV of the 
DM subhalos for both the cold and warm WIMP models. 
With a factor of ^ 20 times larger of the cross section 
for warm WIMP scenario, we can see that the fltixes of 
the most luminous subhalos in the two scenarios are com- 
parable. Also shown in Fig. [7] are the upper limits (for 
80%66-f20%T+T~ case) of dwarf galaxies derived through 
11-month observations of Fermi |2|. The upper limits are 
in general higher than the model expected fluxes, which 
means the first year Fermi data may not be able to probe 
the DM subhalos of both the cold and warm WIMP mod- 
els discussed here. 

Fig. |8] gives the results of the accumulative num- 
ber versus the detection significance, defined as cr = 
-^sig/-\//Vbkg, for i? > 10 GeV and 5-yr exposure of 
Fermi. Here the emission from subhalos within ^haif, 
angular radius containing half of the annihilation lumi- 
nosity, is taken into account. The sky range to calcu- 



^ http://fcrmi.gsfc.nasa.gov/ssc/ 



^ See XTM for a claim of DM signal in the most central region of 
the Galactic center. However, the background and possible point 
source contamination need to be carefully studied. 
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FIG. 7: Integral fluxes above 100 MeV of the DM subhalos 
for the cold and warm WIMP models. The arrows show the 
upper limits of dwarf galaxies given by Fermi observations Q] • 



late the background number of events is adopted to be 
max(0haif,^res), whcre ^rcs ~ 0.1°, is the angular resolu- 
tion of Fermi-LAT at S > 10 GeV fzl. 
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FIG. 8: The accumulative number of subhalos with signifi- 
cance higher than a, for energies _E > 10 GeV and 5-yr expo- 
sure of Fermi-LAT for the cold and warm WIMP scenarios. 



Similar with Fig. 2] the number distribution for warm 
WIMP is flatter than that for cold WIMP. This is a sig- 
nature to distinguish these two scenarios. It is interesting 
to note that the potential detectability for warm WIMP 
might be a little bit better (for high <t ones) than that of 
cold WIMP, although the number of subhalos are signif- 
icantly less. This is because the allowed cross section for 
warm WIMP could be larger in principle. However, it is 
generally difficult to detect the SUSY DM signals from 
subhalos with the Fermi detector, either for the cold or 
the warm WIMP scenarios. 



V. CONCLUSIONS 

Since more and more evidence shows that the DM 
tends to be "warm" instead of "cold" (e.g., [la]), it 
is necessary to investigate the possible consequence on 
DM detections if it is indeed warm. For the canoni- 
cal light WDM particle like the sterile neutrino, most 
of the present DM detection experiments will be useless. 
Alternatively the non-thermally produced warm WIMP 
scenario 1^, [20| might be interesting enough, for both 



the cosmological structure formation and the detection 
of DM particles. The large free-streaming of the DM 
may help to solve the problems of CDM scenario at small 
scale, and the WIMP particles are able to be detected 
with most of the experiments searching for DM. 

Based on the high resolution numerical simulations of 
WDM structure formation, we study the possible 7-ray 
signals from the annihilation of warm WIMPs in the 
Milky Way. The Aquarius CDM simulations are also 
employed to compare with the WDM results. We in- 
vestigate the expected skymaps of the DM annihilation, 
as well as the statistical properties of the subhalos. The 
detectability with Fermi telescope is also discussed for 
two benchmark SUSY models of warm and cold WIMP 
scenarios respectively. Unfortunately we find that the 
detectability of the warm WIMPs with current 7-ray ex- 
periments is very poor. Nevertheless, it is interesting to 
investigate the theoretically expected signatures of the 
7-rays from warm WIMP annihilation, in case that they 
might be detected in future. 

The major conclusions of this work can be summarized 
as follows. 

• Due to a suppression of structure formation in 
WDM scenario, subhalo is much less abundant in 
WDM scenario, resulting in a flatter accumula- 
tive number distribution of J-factor and a differ- 
ent A^(> J) vs. J relation between warm and cold 
WIMP models. 

• We find it is difficult to detect the subhalos with 
Fermi telescope, both for cold and warm WIMP 
scenarios. It is found that the detectablity of 
warm WIMP could in principle be better than cold 
WIMP, because a moderately larger annihilation 
cross section is allowed for warm WIMP scenario, 
with a non-thermal production mechanism [2J] . 

• For indirect WIMP search strategy, the Galactic 
center would likely be prior to dwarf galaxies if 
DM is made of warm WIMPs. For cold WIMPs 
the 7-ray emission due to dark matter annihilation 
from the Galactic center is polluted by the high 
background and the subhalos have been believed to 
be better targets for DM indirect searches. In the 
warm WIMP case, however, the emission from the 
Galactic center could be enhanced due to a larger 
cross section, while the emission from dwarf galax- 
ies is not as significantly enhanced because of the 



decrease of the central DM density and concentra- 
tion. For our benchmark models, the signal of the 
warm WIMP annihilation from the Galactic cen- 
ter will be ~ 20 times stronger than that of cold 
WIMPs, while it is comparable for subhalos. This 
might lead to a different detection strategy in case 
that WIMP is warm. 
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Appendix A: Relic density of non-thermal DM from 
cosmic string decay 

Here wc briefly discuss the non-thermal DM density 
from cosmic string decay. We assume the correlation 
length scale of the string is ^{t) in the friction domi- 
nated epoch. It can be given by ^{t) = ^{tc){t/tc)^/'^ 
[74| . where the initial length <^(<c) ~ -^^^'7^^, A is the 
scalar self-quartic coupling. The production of cosmic 
string loops induce the energy lose of long strings. The 
number density of loops created by long strings can be 
evaluated by [7^, UM 



we find tp ^ X 



Then the reduced 



number density of non-thermal DM particles from decays 
of cosmic string loops can be derived as [l^, ^M 



NT ^ ^^^;,|r- 

A /7T- 



9t~ 9t^9tf 






(A3) 



where Y^ is defined as Y^ 



An 



< 



-4de 
dt' 



(Al) 



njs, s = 2TT^g^T^/A5 is 
the entropy density, g* is effective degrees of freedom at 
the corresponding time. The DM relic density is related 
to Y by rih^ = 2.82 x 10^Y^{m^/GcY). If the DM is 
dominated by non-thermal production, we can get the 
corrected DM relic density fth^ ^0.11 easily by choosing 
I/, A - 1, r - 10^* ~ 100, m^ - 0(10^) GeV and 
Tc-O(109) GeV |li,[2|. 



Appendix B: DM distribution from Aquarius 
simulation 

Here we give the statistical results used for the extrap- 
olation of the unresolved subhalos, based on the resolved 
subhalos from Aquarius CDM simulations. More results 
can be found in Refs. 0, H- 

We bin the luminosities Li of the subhalos with re- 
spect to mass and radius. The left panel of Fig. [5] shows 
the differential distribution of luminosity versus subhalo 
mass, dL/dM, and the right panel shows the spatial dis- 
tribution of the luminosity, dL/dV. When doing this 
analysis we assume that dL/dM is independent with the 
spatial distribution dL/dV @], so that we can use all the 
subhalos to derive both dL/dM and dL/dV. The results 
for WDM are also shown in Fig. |9] for comparison. 

The luminosity-mass relation dL/dM can be fitted 
with a power-law function 



dL/dM oc A/ 



-1.14 



(Bl) 



where i/ is a constant of order 1 
contributes N DM particles. 

Here we only consider the non-thermal DM parti- 
cles from the decay of loops below the temperature T^ 



We assume each loop -y^^ j,^^-^ infer the cumulative luminosity distribution as 



-0.14 



For Mth <C M„ 



we 



(the corresponding time is tchi). For mx ~ 100 GeV, derived in 
T^ ~GeV. Then we can get the DM number density 
by integrating the red-shifted cosmic string loop number 
density 



L(> A/th) « M,-^"-'* - M, 

have L{> Mth) oc M^"^°". Note that this resuh is flatter 
than the mass dependence of the cumulative luminosity 
(oc M^"j^°-22^). This may be due to the 



threshold effect when A/th is close to M^ax ~ 10 Mq. 
For the spatial distribution of the luminosity dL/dV we 
use an iso-thermal l3 function 



NT 



(to) 



Nv 



ia 



atp) 



1) rM. 



(A2) 



dL/dV 



1 



[1 + {r/r,)P] 



(B2) 



where tp is the time at which cosmic string loops which 
are decaying at t^ form. Since the loop density decreases 
sharply with time, we can see the DM density is mainly 
contributed by loops which decay right after i^. It means 
the most of non-thermal DM particles are created at t^ 
instantaneously. 

According to the average radius of loop (formed at 

1 *2 13 

tp) R{tp) ~ X^g^^GfiM^itp, and the loop shrink rate 
dR/dt = -TG^i (F is a constant - 10-20) [zl llil, 



to fit the simulation results. The fitting parameters are 
re « 54 kpc and /3 « 2.76. 

The unresolved subhalos in the CDM simulation is 
derived according to the fitting results of dL/dM and 
dL/dV. The masses of unresolved subhalos are assumed 
to extend to Afmin ~ 10"^ Mq from M^^s « 3 x 10^ Mq. 

The J-factor for unrcsolvc subhalos is 



j:iw 



1 



P%Rq 




(B3) 



10 
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FIG. 9; Differential luminosity-mass relation dL/dM (left) and spatial density of luminosity (right) for subhalos. Red circles 
are for Aq-A-2 CDM simulation, and black triangles are for Aq-AW-2 WDM simulation. Solid lines are the fitting results for 
CDM. 




FIG. 10: Skymaps of the J- factors of the smooth halo (top- left), resolved subhalos (top-right), unresolved subhalos (bottom- left) 
and the total contribution (bottom-right) for CDM. 

Fig. [TU] shows the skymaps of J-factors of the smooth (bottom-right) for CDM. This figure is a reproduction of 
halo (top- left), resolved subhalos (top-right), unre- the result given in Ref. [^. 
solved subhalos (bottom-left) and the total contribution 
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